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Abstract 

The  width  and  energy  of  low-index  interfacial  boundaries  (IFBs)  in  Ni-N^Al  are  calculated  using  first-principles  methods  for  tem¬ 
peratures  ranging  from  0  to  1300  K.  The  low-temperature,  coherent  and  chemically  sharp  (1 00),  (1 1 0)  and  (111)  IFBs  are  studied  using 
conventional  spin-polarized  density  functional  methods.  Cluster  expansion  methods,  as  implemented  in  the  ATAT  software  suite,  are 
used  to  estimate  the  interfacial  excess  free  energies  (lEFEs)  and  composition  and  long-range  order  profiles  of  these  defects  as  a  function 
of  temperature.  The  simple  face-centered  cubic-based  cluster  expansion  produces  interfacial  widths  in  the  range  of  1. 5-3.0  nm  at  1000  K. 
Interfacial  widths  double  in  size  with  an  increase  in  temperature  of  500  K.  We  also  find  that  the  lEEEs  for  the  (100),  (110)  and  (111) 
lEBs  are  strongly  temperature  dependent,  decreasing  by  90%  as  temperature  increases  from  0  to  1000  K.  While  vibrational  and  electronic 
entropic  contributions  were  also  considered,  changes  in  free  energy  are  dominated  by  the  configurational  entropy.  The  predicted 
high-temperature  IEEE  is  approximately  10mJm“^  which  is  in  excellent  agreement  with  previous  fits  to  experimentally  measured 
coarsening  rates. 

Published  by  Elsevier  Ltd.  on  behalf  of  Acta  Materialia  Inc. 
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1.  Introduction 

The  Ni-based  superalloys  have  been  integral  to  the  suc¬ 
cess  of  high-performance  turbine  engines  used  in  aerospace 
propulsion  and  terrestrial  power  generation.  The  high-tem- 
perature  properties  of  these  alloys  reflect  the  contributions 
of  a  combination  of  precipitation  and  solid-solution 
strengthening  mechanisms.  The  properties  of  commercial 
superalloys  are  generally  adjusted  through  a  series  of  heat 
treatments  designed  to  produce  an  optimal  precipitate 
micro  structure.  These  micro  structures  consist  of  a  face- 
centered  cubic  (fee)  Ni  solid-solution  matrix  (y)  and  a 
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distribution  of  ZI2  (7')  precipitates  with  a  volume  fraction 
ranging  from  40%  to  80%  depending  on  the  process  route 
and  intended  application.  The  mechanisms  thought  to 
influence  work  hardening  and  strengthening  are  controlled 
in  part  by  the  nature  of  the  interfacial  boundary  (IFB) 
between  the  solid-solution  Ni  matrix  and  the  ordered-inter- 
metallic  NhAl  precipitate.  The  deformation  (i.e.  yield 
stress,  creep,  fatigue)  of  these  alloys  is  strongly  influenced 
by  the  manner  in  which  dislocations  interact  with  the  pre¬ 
cipitates  and  the  channels  between  them.  In  order  to  model 
the  dislocation-precipitate  evolution  it  is  important  to  have 
a  detailed  understanding  of  the  nature  of  the  y-y'  interface 
over  a  wide  range  of  temperatures.  A  precise  micro  struc¬ 
tural  evolution  method,  one  that  could  predict  the 
size  distribution  and  spatial  correlation  of  precipitates 
during  non-isothermal  processing,  would  be  of  great  value 
to  the  superalloy  community.  Models  of  micro  structure 
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evolution,  such  as  the  phase-field  method,  are  dependent 
on  the  availability  of  reliable  kinetic  as  well  as  bulk  and 
interfacial  thermodynamics  data.  This  includes  the  interfa¬ 
cial  width  and  free  energies  over  the  composition  and 
temperature  range  of  interest. 

The  primary  method  for  estimating  IFB  energies  is 
based  on  the  measurement  of  y'  coarsening  rates.  The  aver¬ 
age  interfacial  energy  over  the  precipitate  surface  is  esti¬ 
mated  by  applying  the  Lifshitz,  Slyozov  and  Wanger 
(LSW)  theory  of  Ostwald  ripening  [1,2]  and  extensions 
thereof.  The  experimental  technique  is  difficult  as  accurate 
measurements  of  average  particle  size  over  time  is  challeng¬ 
ing.  Using  this  technique  early  work  by  various  groups  esti¬ 
mated  the  interfacial  energy  (cr)  as  ranging  from  19  to 
6  mJ  m“^  for  aging  temperatures  ranging  from  773  to 
1073  K,  respectively  [3-6].  In  addition,  the  coarsening 
model  is  not  well  suited  for  concentrated  solid  solutions 
interacting  with  an  intermetallic  precipitate  [7].  Calderon 
et  al.  [7]  extended  LSW  theory  to  incorporate  both  a  mod¬ 
ified  Gibbs-Thomson  equation  to  include  the  effects  of 
non-ideal  solution  thermodynamics  and  the  effects  of  sol¬ 
ute  solubility  in  each  phase  on  the  flux  conservation  condi¬ 
tion  near  the  interface.  The  new  analysis  produced 
corrected  interfacial  energies  ranging  from  42  to  80  mJ/ 
m“^  for  temperatures  ranging  from  848  to  863  K.  This 
advance  also  led  Ardell  to  reassess  existing  coarsening  mea¬ 
surements  in  the  mid-1990s,  producing  a  values  ranging 
from  8.7  to  0.9  mJ  m“^  for  temperatures  from  823  to 
973  K  [8]. 

There  is  significant  scatter  in  the  measured  values  of  a 
and  the  result  is  strongly  dependent  on  what  model  is  used 
to  describe  the  particle  size  distribution  [7].  Further,  recent 
calculations  employing  an  atomistic  scheme  have  demon¬ 
strated  that  at  high  temperatures  the  IFB  spreads  over  a 
wider  region  and  roughens[9].  Ardell  and  Ozolins  [10]  have 
proposed  that  in  this  case  the  rate  limiting  process  for 
coarsening  is  the  diffusion  rate  through  the  diffuse  IFB. 
This  implies  that  coarsening  rates  will  not  have  the  classical 
dependence  but  instead  will  exhibit  a  ^^2  dependence. 
The  new  interpretation  of  coarsening  rates  produces  a  of 
20  mJ  m“^  as  compared  to  29  mJ  m“^  using  the  extended 
LSW  theory  on  the  same  coarsening  data  in  Ni-Al  for 
annealing  temperatures  of  ?^1170  K. 

Independent  first-principles  methods  are  of  interest  for 
this  problem,  as  they  provide  a  framework  for  predicting 
the  relevant  interfacial  widths,  roughness  and  free  energies 
of  the  low-index  IFBs.  By  alloying  these  materials  the  lat¬ 
tice  mismatch  between  Ni  and  Ni3Al  can  be  controlled  and 
the  lattice  mismatch  in  typical  superalloys  is  less  than  0.1%. 
In  turbine  airfoils  the  IFB  between  precipitate  and  the 
matrix  phase  is  often  considered  to  be  coherent.  Such  a 
small  misfit  implies  that  simulation  supercells  without  mis¬ 
fit  dislocations  provide  a  realistic  approximation  of  the  pre¬ 
cipitate-matrix  IFB. 

In  the  earliest  first-principles  study  of  this  class  of  inter¬ 
face  Price  and  Cooper  [11]  used  a  full  potential  linearized 
muffin  tin  orbital  (LMTO)  method  to  estimate  the  (100) 


interfacial  energy  in  the  atomistically  sharp,  zero-tempera¬ 
ture  limit.  This  study  demonstrated  that  there  were  signif¬ 
icant  differences  for  (\^^)Ni/NhAl  interfacial  energies 
treated  in  the  nonmagnetic  and  ferromagnetic  limits.  Treat¬ 
ing  the  interface  using  a  spin-averaged  (nonmagnetic)  local 
density  approximation  (LDA)  functional  method,  they 
found  that  the  interface  had  a  lower  energy  than  the  sur¬ 
rounding  phases  (e.g.  a  negative  interfacial  energy).  This 
study  also  showed  that  spin-polarized  LDA  calculations 
produced  positive  interfacial  energies  in  the  range  50- 
100  mJ  m“^,  thus  spin-polarization  was  necessary  to  pro¬ 
duce  a  stable  interface  between  Ni  and  NhAl.  These  results 
were  independently  confirmed  by  Wolverton  and  Zunger 
[12],  who  found  a  NijAl  phase  that  is  stabilized  relative  to 
Ni  and  Ni3Al  in  the  spin-averaged  approximation.  The 
geometry  of  a  sharp,  stoichiometric  interfacial-boundary 
produces  short-range  order  similar  to  the  My^/  phase  such 
that  the  negative  value  of  the  interfacial  energy  obtained  in 
spin-averaged  calculations  correlates  with  the  predicted 
stability  of  the  NiyAl  phase.  Prior  to  Wolvertons’s  work 
Schonfeld  et  al.  attempted  to  form  this  phase  by  annealing 
a  Ni-rich  Ni-Al  sample  at  775  K  for  24  h  followed  by  a 
rapid  quench.  Subsequent  diffuse  neutron  scattering  mea¬ 
surements  were  unable  to  identify  the  presence  of  the 
NiqAl  phase  [13]. 

In  this  work  we  revisit  the  sharp  interface  problem  using 
a  pseudopotential  plane-wave  method  within  a  spin-polar¬ 
ized  generalized  gradient  approximation  (S-GGA).  In  addi¬ 
tion,  we  explore  the  Ni-rich  part  of  the  binary  Ni-Al  phase 
diagram  using  recently  developed  methods  for  estimating 
the  entropic  contributions  to  the  alloy  free  energy.  These 
methods  are  implemented  in  a  software  package,  the  Alloy 
Theory  Automated  Toolkit,  which  has  been  developed  by 
one  of  the  authors  [14,15].  This  method  is  used  to  evaluate 
the  interfacial  width  and  excess  free  energy  as  a  function  of 
temperature. 

2.  Computational  methods 

A  combination  of  methods,  most  of  which  have  previ¬ 
ously  been  described  in  the  literature,  are  used  to  calculate 
the  structure  and  energies  of  IFBs.  The  Vienna  Ab-initio 
Simulation  Package  (VASP)  is  used  throughout  this  work 
to  calculate  the  structural  properties,  the  formation  energy 
of  various  Ni-Al  ordered  phases  and  the  energies  of  chem¬ 
ically  sharp  interfacial  yjy'  boundaries  [16].  Spin  polariza¬ 
tion,  pseudopotential  and  exchange  correlation  methods 
were  assessed  for  predicting  structural  properties  for  Ni 
and  Ni^Al,  including  lattice  mismatch.  In  general  a  plane- 
wave  cutoff  of  302  eV  was  employed  with  a  k-point  sam¬ 
pling  density  of  8-8-8  for  the  smaller  supercells. 
Numerical  k-space  integration  throughout  the  Brillouin 
zone  is  accomplished  through  a  Gaussian  broadening 
parameter  of  0.2  eV.  Ultrasoft  pseudopotentials  [17]  and 
projected  augmented  wave  (PAW)  [18]  pseudopotentials 
constructed  using  the  local  density  approximation  (LDA) 
[19]  and  the  generalized  gradient  approximation  (GGA) 
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[19]  are  used  to  evaluate  the  lattice  parameters,  lattice  mis¬ 
match  and  elastic  constants  for  Ni  and  Ni3Al  in  order  to 
establish  the  relative  accuracy  of  the  methods.  Calculations 
for  spin-polarized  GGA  employed  the  Vosko,  Wilk  and 
Nusair  interpolation  for  the  correlation  part  of  the 
exchange  correlation  functional  [16,20]. 

Given  the  variety  of  approximations  we  assume  that  the 
most  accurate  representation  of  the  two-phase  system  is 
that  which  produces  the  best  overall  correspondence  to 
experimental  lattice  constants,  elastic  moduli  and  particu¬ 
larly  the  lattice  misfit.  The  results  of  these  calculations 
are  shown  in  Tables  1  and  2.  The  calculated  results  are  sys¬ 
tematically  more  accurate  for  these  materials  using  the 
GGA,  with  a  slightly  better  result  for  spin-polarized  calcu¬ 
lations.  For  calculations  of  the  IFBs,  ultrasoft  pseudopo¬ 
tentials,  within  the  spin-polarized  GGA,  produce  both 
the  lowest  errors  in  the  elastic  constants  and  lattice  misfit. 
This  method  is  then  used  to  calculate  the  chemically  sharp, 
low-temperature  interfacial  y/y'  boundaries  for  three  low- 
index  boundaries:  (100),  (110)  and  (111).  Similarly  these 
approximations  are  used  to  construct  the  reference  config¬ 
urations  for  the  cluster  expansion. 

The  use  of  a  spin-polarized  theory  for  the  intended  pur¬ 
pose  of  modeling  high-temperature  y  and  y'  phases  may  at 
first  seem  counterintuitive.  The  Curie  temperatures  for  Ni 
and  Ni^Al  are  approximately  633  and  48  K  respectively 
[21,22].  However,  the  use  temperatures  for  Ni-based  super¬ 
alloys  can  reach  temperatures  above  1050  K.  Calorimetry 
data  on  the  heat  capacity  of  elemental  Ni,  as  it  passes 
through  the  magnetic  transition  temperature,  indicate  that 
it  likely  retains  large  magnetic  fluctuations.  Further,  results 
of  spin-resolved  X-ray  photoelectron  spectroscopy  suggest 
that  the  spin  polarization  within  0.075  nm  of  the  Ni  core 
does  not  change  with  temperatures  up  to  863  K  (1.37rc) 
[21].  These  observations  suggest  that  as  the  temperature 
increases,  magnetic  contributions  to  the  enthalpy  of  Ni 
remain  important  even  above  the  Curie  temperature.  From 
these  observations  and  previous  calculations  we  conclude 


that  the  ferromagnetic  state  produces  the  best  approxima¬ 
tion  of  the  Ni-NhAl  IFBs  within  the  standard  LDA/GGA 
approximations. 

2.7.  Low-temperature  sharp  interfaeial  boundaries 

Unlike  planar  faults  occurring  in  single-phase  materials 
(i.e.  stacking  faults)  the  interfacial  simulation  cells  include 
strain  energy  produced  by  the  lattice  mismatch  between  the 
two  phases.  Here  we  assume  that  the  lattice  mismatch  is 
small  enough  so  that  the  boundaries  are  coherent,  with 
no  misfit  dislocations  present.  The  strain  energy  needs  to 
be  accounted  for  in  order  to  recover  the  interfacial  energy. 
The  simplest  technique  is  to  calculate  the  interfacial  energy 
for  increasing  cell  wavelengths  (2)  normal  to  the  boundary. 
In  the  low-temperature  limit  the  supercells  consist  of  pure 
Ni  and  stoichiometric  NhAl  blocks,  of  length  2/2,  stacked 
normal  to  the  interface  plane.  The  strain  energy  in  the  cells 
increases  linearly  with  cell  length  and  the  strain-free  inter¬ 
facial  energy  can  be  recovered  by  extrapolating  the  cell 
energy  to  zero  wavelength.  However,  this  technique  is  not 
very  precise,  particularly  when  the  IFB  produces  charge 
perturbations  over  several  layers  normal  to  the  interfacial 
plane.  Previous  studies  have  found  this  approach  unsuit¬ 
able  for  producing  precise  yjy'  IFB  energies  [11]. 

The  second  method  requires  calculating  the  energy  of 
several  strained  reference  cells.  Here  also  the  electronic 
structure,  lattice  positions  and  lattice  parameters  are  opti¬ 
mized  for  interfacial  supercells  of  increasing  wavelength. 
For  an  interfacial  cell  of  length  2  reference  Ni  and  NhAl 
cells  of  length  kjl  are  constructed  that  incorporate  the 
in-plane  lattice  vectors  (e.g.  in  the  plane  of  the  IFB)  derived 
from  the  interfacial  supercell.  Next  the  length  of  these  cells, 
normal  to  the  boundary,  are  optimized  while  holding  the 
in-plane  lattice  vectors  constant.  This  constraint  produces 
an  excellent  approximation  of  the  strain  energy  produced 
in  the  interfacial  cell.  The  interfacial  energy  is  recovered 
by  subtracting  the  energy  of  the  two  references  cells  from 


Table  1 

Calculated  and  measured  [37-39]  lattice  parameters,  elastic  constants  (Mbar)  and  lattice  misfit  (^)  for  Ni  and  NisAl  calculated  using  the  PAW  method. 


PAW 

Metal 

Property 

Spin  averaged 

Spin  polarized 

Exp. 

LDA 

GGA 

LSDA 

S-GGA 

Fee  Ni 

ao  (nm) 

0.34197 

0.35150 

0.34258 

0.35219 

0.35238 

K 

1.144 

1.968 

1.175 

1.942 

1.860 

Cii 

3.185 

2.480 

3.476 

2.704 

2.481 

C12 

2.227 

1.712 

2.019 

1.561 

1.549 

C44 

1.383 

1.121 

1.618 

1.294 

1.242 

(err) 

22% 

5% 

26% 

4% 

LI2  NisAl 

ao  (nm) 

0.34823 

0.35685 

0.349276 

0.35699 

0.35670 

K 

2.183 

1.779 

2.171 

1.773 

1.739 

Cii 

2.774 

2.264 

2.787 

2.343 

2.245 

C12 

1.888 

1.537 

1.863 

1.488 

1.486 

C44 

1.441 

1.189 

1.488 

1.248 

1.244 

(err) 

18% 

0.5% 

18% 

1.4% 

0.01814 

0.01511 

0.01936 

0.01354 

0.01218 

error 

49% 

24% 

59  % 

11% 
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Table  2 

Calculated  and  measured  [37-39]  lattice  parameters,  elastic  constants  (Mbar)  and  lattice  misfit  (^)  for  Ni  and  N^Al  calculated  using  ultrasoft 
pseudopotentials. 


USPP 

Metal 

Property 

Spin  averaged 

Spin  polarized 

Exp. 

LDA 

GGA 

LSDA 

S-GGA 

Fee  Ni 

ao  (nm) 

0.34294 

0.35258 

0.34221 

0.35337 

0.35238 

K 

2.515 

1.985 

2.383 

1.962 

1.860 

Cii 

3.154 

2.506 

3.034 

2.380 

2.481 

C12 

2.195 

1.725 

2.057 

1.753 

1.549 

C44 

1.358 

1.057 

1.363 

1.253 

1.242 

(err) 

20% 

6.8% 

17% 

4.6% 

LI2  NfiAl 

ao  (nm) 

0.34893 

0.35769 

0.34928 

0.35784 

0.35670 

K 

2.163 

1.777 

2.159 

1.787 

1.739 

Cii 

2.749 

2.271 

2.778 

2.397 

2.245 

C12 

1.870 

1.529 

1.849 

1.481 

1.486 

C44 

1.431 

1.173 

1.471 

1.240 

1.244 

(err) 

17% 

0.2% 

18% 

1.8% 

0.01731 

0.01438 

0.02044 

0.01257 

0.01218 

error 

42% 

18% 

68% 

3% 

the  energy  of  the  interfacial  cell  and  normalizing  to  a  stan¬ 
dard  unit  area. 

The  lattice  vectors  used  to  construct  (100),  (110)  and 
(111)  IFBs  are  presented  in  Table  3.  Care  was  taken  in 
maintaining  a  constant  k-point  sampling  density  in  recipro¬ 
cal  space  as  the  aspect  ratios  of  the  cells  change  with 
increasing  wavelength.  A  similar  strategy  was  employed 
for  the  Ni  and  N^Al  cell  reference  calculations  described 
in  the  previous  paragraph. 

2.2.  Cluster  expansion  method  and  the  Alloy  Theoretic 
Automated  Toolkit 

We  rely  on  the  cluster  expansion  formalism  [23-26], 
which  represents  the  energy  of  a  crystalline  alloy  with  a 
computationally  efficient  Hamiltonian  taking  the  form  of 
a  polynomial  in  terms  of  occupation  variables  ct/  =  ±1 
indicating  the  type  of  atom  (Al  or  Ni)  residing  on  each  lat¬ 
tice  site  i\ 

E  =  yAijOiOj  +  y^JijkOiajak  H - 

The  unknown  coefficients,  J  .,  of  this  polynomial  are  called 
effective  cluster  interactions  (ECIs)  and  are  fit  to  a  database 
of  ab  initio  structural  energies.  It  has  been  formally  shown 
that  an  infinite  series  of  this  form  can  represent  any 

Table  3 

Supercell  geometries  used  to  calculate  the  low-temperature  sharp  interfa¬ 
cial  energies.  Lattice  vectors  are  in  units  of  the  lattice  parameter  and  the 
scaling  parameter  (N)  is  varied  in  order  to  sample  different  cell 
wavelengths  normal  to  the  IFB.  The  k-point  sampling  is  reduced  with 
increasing  wavelength  in  order  to  maintain  a  constant  k-space  sampling 
density. 


Boundary 

Lattice  vectors 

a  b 

C 

Scaling  factor  (N) 

(001) 

(100) 

(010) 

(OOl)N 

2,  4,  8,  16 

(110) 

(1 10) 

(001) 

(IIO)N 

2,  4,  8,  16 

(111) 

(1 10) 

(112) 

(lll)N 

2,  4,  8,  16 

configuration  dependence  of  the  energy  [23].  Moreover,  for¬ 
mal  methods  have  been  developed  to  determine  the  number 
of  terms  and  the  database  size  needed  to  achieve  a  given  pre¬ 
cision  [14].  The  ECIs  can  typically  be  determined  from  a  rea¬ 
sonably  sized  database  of  ab  initio  calculations  involving 
small-unit-cell  atomic  arrangements.  These  ECIs  can  then 
be  used  in  a  large-scale  equilibrium  lattice-gas  Monte  Carlo 
simulations  of  the  coherent  interface,  without  necessitating 
repeated  large-scale  ab  initio  calculations  for  each  atomic 
configuration  visited  in  a  thermodynamic  equilibrium.  These 
tasks  were  carried  out  with  the  help  of  the  Alloy  Theoretic 
Automated  Toolkit  (AT AT)  [14,15,27,28]. 

3.  Results 

3.1.  Low-temperature  sharp  interfacial  boundaries 

Optimizing  the  electronic  structure  of  high-aspect-ratio 
supercells  can  be  challenging.  Avoiding  charge  oscillations 
in  the  long- wavelength  direction  during  band-by-band  opti¬ 
mization  can  be  difficult  even  using  the  iterative  diagonaliza- 
tion  methods.  Results  representative  of  the  low- temperature 
sharp-interface  limit  for  (100),  (1 10)  and  (1 1 1)  y-y'  IFBs 
(i.e.  pure  Ni  and  stoichiometric  N^Al)  are  shown  in  Figs.  1 
and  2.  The  interfacial  energies  are  very  sensitive  to  errors 
in  the  local  spin  moments.  In  order  to  illustrate  the  conver¬ 
gence  of  the  electronic  structure  as  a  function  of  cell  wave¬ 
length  (2)  the  total  magnetic  moment  of  the  interfacial  and 
Ni  and  NLAl  reference  cells  are  shown  in  Fig.  1 .  As  expected 
the  supercell  moment  increases  linearly  with  increasing  2  and 
the  increase  in  moment  with  cell  size  is  the  same  for  the  defect 
and  reference  cells.  Also  the  cells  containing  the  IFB  have  a 
slightly  smaller  moment  than  the  reference  cells  reflecting  a 
slight  decrease  in  the  Ni  moment  near  the  IFBs.  This  obser¬ 
vation  is  consistent  with  that  reported  by  Price  and  Cooper 
using  the  full  potential  LMTO  method  [11]. 

Results  for  the  sharp  interfacial  energies  produced  using 
strained  Ni  and  N^Al  reference  cells  are  shown  in  Fig.  2. 
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Wavelength  (Lattice  Units) 


Fig.  1 .  Magnetic  moments  as  a  function  of  wavelength  for  supercells  used 
to  calculate  the  interfacial  energies.  Filled  markers  and  dashed  lines 
correspond  to  the  supercells  containing  the  IFBs  and  open  markers  and 
solid  lines  correspond  to  the  reference  cells  for  pure  Ni  and  stoichiometric 
Ni^Al.  The  latter  magnetic  moments  are  calculated  by  taking  the  average 
of  Ni  and  Ni^Al  cells  the  same  wavelength  as  the  defect  cell.  As  expected 
the  moments  scale  with  the  length  of  the  cells,  and  the  interfacial  and 
reference  cells  have  nearly  identical  moments.  The  units  for  the  length  of 
the  cells  are  taken  as  the  periodic  length  in  the  (100),  (110)  and  (111) 
directions:  ao,  V^slq  and  V^slq,  respectively,  in  lattice  units. 


Wavelength  (Lattice  units) 

Fig.  2.  The  zero-temperature  sharp  interfacial  energies  as  a  function  of 
wavelength  for  the  (100),  (110)  and  (111)  Ni-Ni3Al  IFBs  (i.e.  pure  Ni 
and  stoichiometric  Ni^Al).  Previous  full  potential  LMTO  result  of  Price 
and  Cooper  for  the  (100)  boundary  are  shown  as  open  triangles  [11].  The 
previous  work  did  not  include  optimization  of  the  atomic  positions,  which 
may  explain  the  discrepancy  with  the  current  results.  Convergence  with 
cell  size  is  more  rapid  along  (100)  and  (110)  as  compared  to  (111). 

Included  in  this  figure  are  the  previous  results  of  Price  and 
Cooper  [11]  for  the  (100)  IFB.  The  current  calculations 
include  more  degrees  of  freedom  that  are  optimized  during 
the  relaxation  process,  so  it  is  not  surprising  that  the  cur¬ 
rent  (100)  interfacial  energies  converge  to  a  different  value. 
The  energies  for  the  sharp  (001),  (1 10)  and  (1 1 1)  IFBs 
converge  to  approximately  41,  55  and  48  mJ  m“^  respec¬ 
tively  in  the  long-cell-wavelength  limit. 

3.2.  Temperature-dependent  results  from  Monte  Carlo 
simulations 

A  cluster  expansion  (CEl)  was  developed  for  binary  fee 
structures  in  the  M-rich  portion  of  the  phase  diagram.  This 
was  constructed  from  25  superstructures  between  the  Ni  fee 


and  LI2  endpoints,  including  18  structures  automatically 
generated  by  AT  AT  and  seven  additional  superstructures 
(“special  quasirandom  structures”  [29]  and  dilute  impuri¬ 
ties).  All  the  reference  calculations  were  run  using  ultrasoft 
pseudopotentials  within  the  spin-polarized  GGA  approxi¬ 
mation.  The  ATAT  process  produces  a  cluster  expansion 
with  seven  pairs  within  0.67  nm  and  two  triplets  within 

0.36  nm  with  an  overall  cross-validation  score  of 
0.008  eV.  Here  both  configurational  and  vibrational  entro- 
pic  terms  were  assessed  using  techniques  described  previ¬ 
ously.  The  resulting  effective  cluster  interactions  were 
then  used  to  generate  the  Ni-rich  part  of  the  Al-Ni  phase 
diagram  which  is  shown  in  Fig.  3.  While  the  y'  side  of  the 
7  +  7'  two-phase  region  boundary  is  in  good  agreement 
with  the  experimental  results  of  Zhang  et  al.  (shown  as  a 
solid  line)  [30],  the  predicted  7'  +  7-7  boundary  is  higher 
than  the  experimental  fit.  At  15  at.%  Al  the  predicted  sol- 
vus  line  is  300  K  higher  than  that  interpolated  from 
experiment. 

It  is  worth  noting  that  in  a  preliminary  study  effective 
cluster  interactions  were  developed  using  the  spin-averaged 
GGA  approximation.  The  resulting  phase  diagram  included 
the  additional  NiqAl  phase  anticipated  by  previous  calcula¬ 
tions  based  on  spin-averaged  LDA  [12].  The  phase  lies  in 
the  middle  of  the  composition  range  between  7  and  7'  and 
persists  to  temperatures  above  600  K.  As  mentioned  previ¬ 
ously,  to  our  knowledge,  there  is  no  experimental  evidence 
that  the  NiqAl  phase  is  stable  in  the  binary  Ni-Al  system. 

3.3.  Temperature-dependent  interfaeial  boundary  widths 

The  IFBs  are  simulated  using  high-aspect-ratio  cells, 
where  the  periodic  length  normal  to  the  IFB  is  at  least  twice 


Atomic  Fraction  Ni 

Fig.  3.  The  Ni-rich  phase  diagram  produced  by  the  first  cluster  expansion 
(CEl).  Solid  lines  are  from  the  CALPHAD  calculations  using  the 
PANDAT  method  [30,40]  and  the  dashed  lines  are  cluster  expansion 
predictions  using  Monte  Carlo  methods.  For  CEl  the  ground  state  search, 
and  the  formation  of  the  effective  interaction  coefficients,  was  initialized 
with  the  fee  endpoints  Ni  and  LI2  N^Al.  The  (7  +  /)/7  phase  boundary  is 
shifted  to  higher  Ni  compositions  relative  to  experimental  observations. 
At  Ni-15A1  (at.%),  represented  in  the  figure  as  a  vertical  dotted  line,  the 
difference  between  the  predicted  and  measured  phase  boundary  is 
approximately  300  K. 
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the  interface-plane  periodic  length.  For  example  the  (100) 
IFB  was  approximated  with  a  cell  with  cubic  lattice  vectors 
of  length  32  x  32  x  128.  The  initial  cell  contained  Ni  and 
Ni2,Al  slabs,  of  equal  width,  with  a  sharp  (100)  boundary 
aligned  normal  to  the  longest  lattice  vector  (e.g  the  [001] 
direction).  The  cell  was  equilibrated  at  1000  K  using  a 
two-step  process.  First,  the  cell  composition  had  to  be 
adjusted  to  reflect  our  target  cell-composition  of  15  at.% 
Al.  This  was  accomplished  by  introducing  a  random  distri¬ 
bution  of  Ni  atoms  on  the  Al  sublattice  in  Ni^Al.  Next,  this 
cell  was  equilibrated  using  a  canonical  ensemble  Monte 
Carlo  approach  at  this  composition.  After  reaching  equi¬ 
librium  the  Monte  Carlo  simulation  is  continued  and  snap¬ 
shots  of  the  IFB  were  collected  at  time  intervals  greater 
than  200  steps,  ensuring  uncorrelated  measurements.  For 
each  measurement  the  concentration  and  long-range  order 
(LRO)  parameter  are  averaged  over  planes  normal  to  the 
IFB.  These  profiles  are  then  fit  to  a  hyperbolic  tangent, 
the  individual  profiles  are  shifted  so  that  they  share  a  com¬ 
mon  point  of  inflection  and  then  all  the  data  is  binned  nor¬ 
mal  to  the  IFB.  The  resulting  curve  is  again  fit  to  a  general 
hyperbolic  tangent  function  and  the  10-90  width  is  taken 
from  this  fit.  An  example  of  the  final  (100)  interfacial  pro¬ 
file  obtained  by  this  procedure  is  shown  in  Fig.  4.  Error 
bars  on  the  data  points  represent  one  standard  deviation 
and  are  obtained  from  averaging  over  40  separate  Monte 
Carlo  snapshots.  At  1000  K  the  interfacial  10-90  width  is 
2.06  zb  0.03  and  1.80  zb  0.03  nm  for  the  composition  and 
LRO  parameter,  respectively. 

The  (110)  and  (111)  IFB  widths  were  calculated  using  a 
similar  procedure.  The  temperature  dependence  for  the 
concentration  and  LRO  parameter  profiles  of  the  three 
low-index  boundaries  are  plotted  in  Figs.  5  and  6  respec¬ 
tively.  At  any  given  temperature  the  width  of  the  (111) 
IFB  is  approximately  70%  larger  than  the  other  bound¬ 
aries.  The  rate  of  increase  in  interfacial  width  is  approxi¬ 
mately  0.21  nm  for  every  100  K  for  the  (100)  and  (1 10) 
boundaries,  with  concentration  and  LRO  having  the  same 
trends.  The  rate  of  change  for  the  (111)  boundary  is 
slightly  more  rapid,  0.23  and  0.26  nm  per  100  K  for  con¬ 
centration  and  LRO,  respectively.  For  a  fixed  composition, 
as  the  temperature  approaches  the  solvus  boundary,  the 
Ni^Al  slab  becomes  thinner,  making  it  more  difficult  to  fit 
the  boundary  profiles.  This  is  reflected  in  the  poorer  statis¬ 
tics  found  for  the  higher-temperature  boundaries. 

3.4.  Interfacial  excess  free  energy 

One  of  the  primary  goals  of  this  study  is  to  estimate  the 
y-y'  interfacial  excess  free  energy  (lEFE).  In  order  to 
accomplish  this  the  following  approach  was  used.  For  each 
low-index  IFB  two  cells  were  constructed,  with  one  cell 
having  twice  the  length  (2)  of  the  other  normal  to  the 
respective  IFBs.  These  cells  were  then  equilibrated  over  a 
range  of  temperatures,  typically  from  10  to  1300  K,  for 
an  overall  target  composition  of  15  at.%  Al.  At  each 
temperature  the  equilibrated  cells  have  the  same  volume 


Fig.  4.  The  LRO  and  composition  profile  for  the  (1 00)  y-y'  IFB  at  1000  K 
derived  from  Monte  Carlo  simulations  using  the  first  cluster  expansion 
(CEl).  Each  data  point  in  the  profile  is  derived  from  40  uncorrelated 
Monte  Carlo  snapshots  of  the  IFB  based  on  an  overall  cell  composition  of 
Ni-15A1  (at.%).  The  profiles  are  centered  about  a  common  point  of 
inflection,  the  data  is  then  binned  and  averaged  for  each  plane  normal  to 
the  IFB.  Each  of  the  data  points  represent  one  of  the  (100)  lattice  planes, 
and  includes  error  bars  for  the  expected  numerical  error  in  the  compo¬ 
sition  (i.e.  one  standard  deviation). 
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Fig.  5.  The  low-index,  y-y'  interfacial  10-90  widths,  based  on  the 
composition  profiles,  as  a  function  of  temperature.  The  results  are  based 
on  an  overall  cell  composition  of  Ni-15A1  (at.%)  and  CEl. 


fraction  of  y  and  y'  and  the  interfacial  excess  energy 
(IFF)  is  derived  by  subtracting  the  overall  energy  of  the 
two  cells  {2Ex/2  —  Ex).  Fig.  7  shows  the  lEFE  as  a  function 
of  temperature  for  the  (100)  IFB  as  calculated  using  cluster 
expansion  CEl.  The  low- temperature  boundary  energies 
were  calculated  directly  using  a  conventional  sharp-inter- 
face  model  in  the  previous  section.  Comparing  the  pre¬ 
dicted  energies  from  these  two  techniques  gives  some 
indication  of  self-consistency  between  the  cluster  expansion 
and  the  conventional  sharp-interface  model.  The  cluster 
expansion,  CEl,  produces  energies  of  97,  94  and  40  mJ  m“^ 
for  the  (100),  (110)  and  (111)  IFBs  respectfully.  The 
values  for  the  (100)  and  (110)  IFBs  are  approximately 
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Temperature  (K) 

Fig.  6.  The  low-index  y-y'  interfacial  10-90  widths,  based  on  the  LRO 
parameter,  as  a  function  of  temperature.  The  results  are  based  on  an 
overall  cell  composition  of  Ni-15A1  (at.%)  and  CEl. 

twice  that  found  by  directly  calculating  this  energy  using 
VASP  (see  Section  3.1)  while  the  (111)  energy  is  within 
10  mJ  m~^.  While  these  are  substantial  differences  in  ener¬ 
gies  at  the  low-temperature  limit,  it  remains  to  be  seen  if 
the  lEFE  is  strongly  affected  by  this  at  higher  temperatures. 
In  order  to  assess  this  we  derive  the  change  in  the  IEEE 
with  temperature  and  return  to  this  issue  by  considering 
a  second  cluster  expansion  in  the  discussion  section. 

The  IEEE  is  calculated  by  integrating  the  lEE  over  tem¬ 
perature  using  a  generalization  of  the  Gibbs-Helmholtz 
relation  {dF /T /dT)  =  —E/T^,  where  F  and  E  are  the 
Helmholtz  and  internal  energies,  respectively): 

IEFE{T)  =E^^-T  [  (1) 

Jo  T 

The  IEE{T)  is  approximated  by  fitting  the  interfacial 
excess  energy  as  a  function  of  temperature  at  low  tempera¬ 
ture  (<  300  K)  to  and  at  high  temperatures 

(>  300  K)  to  a  fifth-  or  sixth-order  polynomial  in  T.  The 
resulting  piecewise  continuous  function  is  numerically  inte¬ 
grated  using  Eq.  (1),  yielding  the  IEEE  as  a  function  of 
temperature.  An  example  of  the  excess  energy  variation 


Temperature  (K) 

Fig.  7.  The  lEE  as  a  function  of  temperature  for  the  (100)  y-y'  boundary 
calculated  using  CEl.  Open  circles  are  the  predicted  energies,  and  the  solid 
line  is  a  piecewise  continuous  fit  to  the  data. 


with  temperature  and  the  fit  to  this  function  is  shown  in 
Fig.  7.  The  solid  line  is  an  example  of  a  typical  fit  to  the 
derived  energies. 

Fig.  8  shows  the  IEEE  (i.e.  Eq.  (1)),  derived  using  CEl, 
for  the  low-index  boundaries  at  a  nominal  overall  cell  com¬ 
position  of  15  at.%  Al.  Note  the  rapid  decrease  in  IEEE 
with  increasing  temperature.  Comparing  the  predicted 
energies  to  that  estimated  from  fits  to  coarsening  experi¬ 
ments  [8,10]  we  find  very  good  agreement  with  experiment 
for  the  lEFEs.  The  open  circles  in  the  figure  are  from  fits  to 
the  coarsening  data  assuming  that  the  free  energy  of  mixing 
is  ideal  this  reduces  the  interfacial  free  energy  within  the 
LSW  model  [8].  The  two  filled  markers  in  the  figure  are 
results  from  the  most  recent  fits  to  coarsening  data.  One 
of  the  data  points  was  derived  using  the  trans-interface-dif- 
fusion-controlled  theory  of  coarsening  which  assumes  that 
the  rate-limiting  process  for  coarsening  is  the  diffusion  rate 
through  the  diffuse  (chemistry  and  LRO)  IFB.  Thus  the 
filled  circle  and  square  are  the  result  of  fits  to  the  same 
coarsening  data  using  the  classical  LSW  dependence 
and  a  dependence,  respectively  [10]. 

4.  Discussion 

4.1.  Idealized  low-temperature  boundaries 

The  low-temperature  sharp  IFBs  have  been  previously 
modeled  using  density  functional  theory  by  several  groups. 
Price  and  Cooper  [11]  used  the  full  potential  LMTO 
method  to  explore  the  role  of  spin  polarization  on  the  sta¬ 
bility  of  the  (100)  y-y'  IFB.  Using  spin-averaged  calcula¬ 
tions  on  the  (100)  boundary  yielded  negative  energies 
ranging  from  -23  to  -1 18  mJ  m“^  for  2  =  1  and  3,  a  result 
not  consistent  with  experimental  observations.  We  have 
run  similar  calculations  using  GGA,  with  spin-averaged 
ultrasoft  pseudopotentials,  and  find  for  a  (100)  IFB,  with 


Fig.  8.  The  lEFE  as  a  function  of  temperature  for  the  (100),  (110)  and 
(111)  y-y'  boundaries  calculated  using  CEl.  The  lEFE  is  derived  by 
numerically  integrating  Eq.  (1)  using  the  piecewise  continuous  polynomial 
fit  to  the  interfacial  excess  energy.  The  solid  and  dashed  lines  are  the 
predicted  lEFEs  which  are  in  good  agreement  with  energies  derived  from 
earlier  coarsening  experiments  of  Ardell.  Open  and  filled  circles  are  from 
coarsening  data  fit  to  LSW  theory  [8]  and  the  filled  square  is  from  fits 
using  trans-interface-diffusion-controlled  theory[10]. 
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X  =  12,  an  energy  of  -46  mJ  m“  .  Proceeding  with  spin- 
polarized  calculations  on  the  same  boundary,  Price  and 
Cooper  found  that  the  spin  moments  transitioned 
smoothly  over  2-3  atomic  layers,  and  identified  this  spin 
density  transition  across  the  boundary  as  the  main  contri¬ 
bution  to  the  interfacial  energy.  The  spin-polarized  calcula¬ 
tions  yielded  positive  interfacial  energies  of  mJ  m“^ 
(see  Fig.  2).  We  find  a  similar  transition  of  spin  density 
across  all  three  low-index  IFBs,  and  interfacial  energies 
generally  lower  than  the  initial  calculations  of  Price  and 
Cooper.  This  is  likely  due  to  relaxing  additional  degrees 
of  freedom  over  larger  supercells  (i.e.  atomic  positions 
and  lattice  parameters). 

The  low- temperature  chemically  sharp  IFB  energies  for 
the  (100),  (110)  and  (111)  are  41,  55  and  48  mJ  m“^ 
respectfully  using  density  functional  methods  (ultrasoft 
pseudopotentials  within  S-GGA).  Following  Price  and 
Cooper,  the  relative  size  of  the  interfacial  energies  is  likely 
related  to  the  degree  to  which  the  magnetic  moment  of  Ni 
is  disrupted  by  the  introduction  of  the  IFB  along  the  three 
low-index  directions.  To  a  first-order  approximation,  the 
relative  size  of  the  interfacial  energies  reflects  the  number 
of  atomic  spin  moments  suppressed  within  a  lattice  con¬ 
stant  on  the  Ni  side  the  boundary. 

More  recently  there  has  been  one  other  first-principles 
study  of  the  IFBs  which  focussed  on  tungsten  segregation 
to  the  (100)  y-y'  interface,  Amouyal  et  al.  [31]  used  the 
local  density  approximation  as  implemented  in  VASP  to 
calculate  the  low-index  IFB  energies.  They  report  energies 
of  the  pure  Ni/NhAl  interfaces  of  42.6,  29.1  and 
16.0  mJm“^  for  the  (100),  (110)  and  (111)  interfaces, 
respectively.  These  are  significantly  lower  than  the  con¬ 
verged  values  found  in  the  current  work. 

4.2.  Interfacial  boundaries  at  high  temperatures 

Using  the  cluster  expansion  approach  we  find  the  IFBs 
at  finite  temperatures,  as  described  through  the  composi¬ 
tion  and  LRO  parameter,  change  continuously  across  the 
boundary.  As  temperature  increases  the  boundary  broad¬ 
ens  significantly.  One  interesting  aspect  of  the  interfacial 
widths  is  the  consistently  larger  (111)  IFB  relative  to  the 
other  low-index  boundaries.  First,  this  result  is  consistent 
with  the  low-temperature  sharp  interface  results  in  that 
the  (111)  boundary  converges  at  a  significantly  larger 
wavelength  than  the  other  two  interfaces.  Also,  as  pointed 
out  by  Ardell  and  Ozolins  [9],  the  diffusion  rates  are 
expected  to  be  1-3  orders  of  magnitude  times  slower  in 
the  disordered  NhAl  boundary  than  in  disordered  Ni-Al. 
Currently,  the  cuboidal  structure  of  y'  precipitates  often 
observed  in  the  superalloys  is  attributed  to  the  strain  field 
produced  by  the  lattice  mismatch  between  the  two  phases. 
Combined  with  these  observations,  the  current  results 
imply  that  there  may  be  a  window  of  precipitate  size, 
between  the  critical  size  and  when  the  strain  field  domi¬ 
nates,  where  the  kinetics  would  strongly  favor  precipitate 
growth  along  (100)  and  (110)  as  compared  to  (111). 


The  predictions  for  interfacial  widths  may  prove  useful 
in  models  of  micro  structure  evolution  that  explicitly  treat 
the  transition  between  two  domains  (e.g.  different  phases, 
crystal  orientation)  using  a  regular  grid  of  sampling  points. 
Given  the  predicted  changes  in  interfacial  widths  with  tem¬ 
perature  it  is  likely  that  incorporating  these  effects  will  be 
important  when  modeling  non-isothermal  processes.  For 
Ni-based  superalloys,  solution  heat-treatments  and  in-ser¬ 
vice  temperatures  can  range  from  900  to  1625  K.  This  cor¬ 
responds  to  a  change  in  predicted  boundary  width  of 

1.2  nm,  in  some  cases  doubling  the  width  of  the  boundary. 
In  order  for  models  using  a  regular  grid  to  satisfy  Nyquist’s 
theorem  at  least  two  grid  points  should  fall  within  the 
boundary.  This  requirement  effectively  limits  the  scale  of 
such  calculations  to  the  processing  limitations  of  available 
computer  hardware.  For  example,  in  the  low- temperature 
limit  grid  spacing  should  be  less  than  0.5  nm  to  model 
the  (110)  and  (111)  concentration  profile  in  y-y' .  Assum¬ 
ing  a  practical  upper  limit  on  the  fast  Fourier  transform 
(FFT)  grid  of  512^,  such  methods  applied  to  this  material 
are  limited  to  volumes  of  less  than  a  cube  with  edge  length 
of  0.25  jim  at  900  K.  It  is  not  unusual  to  find  y'  particles 
larger  than  1  pm^  in  engineering  alloys  and  this  could  limit 
the  use  of  FFT-based  micro  structural  evolution  codes  in 
this  class  of  problems. 

Recent  advances  in  electron  optics  have  enabled  a  much 
more  detailed  characterization  of  the  IFB.  There  have  been 
several  recent  studies  employing  high-resolution  scanning 
transmission  electron  microscopy  (HRSTEM)  with 
atomic-resolution  Z  contrast  imaging  (e.g.  high-angle 
annular  dark-field;  HAADF).  Srinivasan  et  al.  [32]  charac¬ 
terized  a  (100)  y-y'  IFB  in  the  Ni-based  superalloy  Rene 
88DT  that  had  been  annealed  at  1039  K  for  50  h.  By  align¬ 
ing  the  sample  along  (0 1 0)  columns  of  atoms  can  be  iden¬ 
tified,  using  image  intensity  and  contrast,  that  correspond 
to  the  y  and  y'  phases.  Using  just  this  information  the  posi¬ 
tion  and  width  of  the  boundaries  was  estimated  as  2.7  and 

1.3  nm  for  the  concentration  and  LRO  parameter,  respec¬ 
tively.  In  order  to  assess  the  compatibility  of  this  approach 
with  the  techniques  used  in  this  paper  we  have  applied  a 
similar  intensity  and  contrast  method  based  on  the  concen¬ 
tration  of  Ni  along  atomic  columns  adjacent  to  a  (100) 
interface,  oriented  as  described  above.  We  find  that  the 
two  methods  produce  similar  IFB  widths  that,  for  the  cases 
tested,  were  within  statistical  error.  While  the  predicted 
interfacial  widths  are  within  the  range  measured  in  the 
superalloy  Rene  88DT,  the  current  calculations  generally 
predict  the  LRO  and  composition  interfacial  widths  to  be 
approximately  equal  over  a  wide  range  in  temperature. 

The  HRSTLM-HAADF  characterization  also  showed 
that  the  midpoint  (i.e.  point  of  inflection)  of  the  composi¬ 
tion  and  LRO  profiles  are  offset  by  approximately 
0.6  nm,  with  the  LRO  midpoint  closer  to  the  y  side  of  the 
boundary  [32].  This  is  a  particularly  interesting  observa¬ 
tion,  because  it  allows  the  system  to  stabilize  shells  of  mate¬ 
rial  around  y'  particles  that  do  not  share  the  ordering  or 
chemistry  of  either  of  the  equilibrium  phases.  Current 
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Monte  Carlo  simulations  of  the  (100)  and  (110)  bound¬ 
aries  based  on  the  cluster  expansion  (CEl)  exhibit  a  small 
shift,  Sc-LRO  =  {arcTanh{0)(.)  —  {arcTanh{^)  ,  which 
reaches  a  maximum  of  — 0.1  nm  at  TOOK.  However, 
the  (111)  boundary  produces  a  positive  shift  that  increases 
to  almost  0.4  nm  at  higher  temperatures,  though  the  error 
in  this  value  also  increases  with  temperature  as  shown  in 
Fig.  9.  Hence,  even  using  a  relatively  simple  binary  cluster 
expansion  to  compute  via  the  Monte  Carlo  method  the 
entropic  contributions  to  the  interfacial  width  can  repre¬ 
sent  the  variation  in  behavior  found  in  complex  alloys  such 
as  Rene  88DT. 

A  second  more  recent  study  by  Ge  et  al.  [33]  performed 
similar  HAADF  imaging  of  a  ( 1 0  0)  boundary  in  a  second- 
generation  Ni-based  superalloy  (DD6).  This  particular 
alloy  contains  2  wt.%  Re,  an  element  not  present  in 
Rene-88DT,  which  is  known  to  improve  high-temperature 
creep  performance  in  this  class  of  alloys.  The  DD6  sample 
had  a  final  anneal  at  1 143  K  for  32  h.  The  HAADF  images 
of  an  (100)  IFB  were  analyzed  using  image  intensity  and 
contrast  to  identify  the  concentration  and  FRO  profiles 
across  the  interface.  In  this  case  the  interfacial  widths 
was  found  to  be  approximately  2.2  nm  for  both  concentra¬ 
tion  and  FRO  parameter,  and  in  contrast  to  earlier  work 
on  Rene  88DT  there  is  no  apparent  shift  in  the  point  of 
inflection  of  the  two  profiles.  Perhaps  fortuitously,  at  these 
annealing  temperatures  the  measured  width  is  in  good 
agreement  with  the  predicted  widths  for  the  (100)  and 
(110)  IFBs  forM-AM/. 

More  recently  Plotnikov  et  al.  [34]  have  used  3-D  atom 
probe  tomography  to  characterize  y'  particles  in  binary  Ni- 
A1  alloys.  This  study  focused  on  the  evolution  of  particle 
size  and  IFB  width  of  during  annealing  at  two  tempera¬ 
tures  (823  and  873  K)  for  times  ranging  from  5  min  to 
4096  h.  The  IFB  widths  were  determined  using  a  spline  fit 
to  the  proximity  histogram  of  the  composition  profile. 
The  10-90  width  was  determined  for  profiles  at  two 
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Fig.  9.  The  shift  in  the  point  of  inflection  {Sc-lro)  for  the  concentration 
and  LRO  proflles  of  the  (111)  y-y'  IFB.  Results  are  based  on  an  overall 
cell  composition  of  Ni-15A1  (at.%)  and  CEl. 


compositions,  Ni-12.5A1  and  Ni-13.5A1.  For  the  former 
composition  aged  at  873  K  for  256  h  the  measured  interfa¬ 
cial  width  is  approximately  1.6  nm;  this  in  good  agreement 
with  the  predicted  values  for  Ni-15A1  at  this  temperature. 
It  is  worth  noting  that  capillary  effects  on  the  chemistry  of 
aged  submicron  precipitates  are  likely  to  produce  Ni-rich  y' 
compared  to  the  slabs  used  in  the  current  calculations. 
However,  the  interfacial  width  is  weakly  dependent  on 
the  Al  content  as  shown  in  the  atom  probe  study  [34]. 

4.3.  Interfacial  excess  free  energies 

The  IFFFs  calculated  using  CEl  drop  precipitously  with 
increasing  temperature.  As  shown  in  Fig.  8,  predicted  val¬ 
ues  at  high  temperature  fall  within  the  range  of  the  limited 
measurements  based  on  growth  kinetics.  These  calculations 
have  also  been  performed  with  and  without  the  vibrational 
and  electronic  entropic  contributions  using  techniques 
described  previously  in  the  literature  [14].  Surprisingly, 
we  find  that  the  free  energy  of  the  diffuse  interfaces  is  dom¬ 
inated  by  the  configurational  entropy.  Recently  Mao  et  al. 
[35]  estimated  the  vibrational  entropic  contributions  to  the 
free  energy  of  sharp  (100),  (110)  and  (111)  y-f  interfaces 
in  Ni-Ni3Al  using  the  harmonic  approximation.  They 
found  that  there  are  significant  vibrational  energy  contri¬ 
butions  to  the  free  energy  for  the  sharp,  ordered  IFBs. 
The  discrepancy  between  these  two  results  is  likely  due  to 
the  differences  in  vibrational  entropy  of  the  high-tempera¬ 
ture  boundary  as  represented  in  the  ordered  and  disordered 
limits.  Previously,  van  de  Walle  et  al.  [36]  have  shown  that 
changes  in  vibrational  entropy  between  ordered  and  disor¬ 
dered  NisAl  is  small  over  a  wide  temperature  range.  In  the 
current  study,  as  the  temperature  rises  the  IFB  becomes 
disordered  and  diffuse,  so  changes  in  the  vibrational 
entropy  contributions  will  be  small.  If  the  boundary 
remains  ordered,  as  the  temperature  increases,  changes  in 
the  vibrational  entropy  may  begin  to  dominate  the  entropic 
contributions  to  the  free  energy.  These  results,  taken  with 
the  insights  from  the  previous  work  of  Mao  et  al.  on  sharp 
interfaces,  may  help  explain  trends  in  free  energies  as 
groups  begin  exploring  the  effects  of  more  complex  chemis¬ 
tries  on  the  energies  of  these  interfaces. 

The  initial  cluster  expansion  of  the  Ni-Al  system  pro¬ 
vides  an  adequate  approximation  of  the  Ni-rich  phase 
diagram  and  the  width  of  the  IFBs  at  elevated  tempera¬ 
tures.  However,  this  cluster  expansion  does  not  produce 
energies  consistent  with  the  low-temperature  sharp-inter- 
facial  energies  found  using  explicit  density  functional  cal¬ 
culations  as  described  in  Section  3.1.  In  order  to 
determine  if  this  has  a  significant  effect  on  the  lEFE  or 
its  temperature  dependence,  a  second  cluster  expansion 
was  considered. 

The  second  cluster  expansion  (CE2)  was  developed  by 
adding  to  the  set  of  input  structures  the  three  shortest- 
wavelength  IFBs  calculated  previously  (see  Section  3.1). 
The  intent  was  to  bias  the  expansion  to  maintain  a  more 
precise  representation  of  the  low-temperature  IFB  energy. 
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Within  the  AT  AT  process  this  produced  four  pairs  within 

0.5  nm,  and  a  triplet  and  quartet  within  0.25  nm  with 
a  cross-validation  score  of  0.008  eV.  The  resulting  effective 
cluster  interactions  were  then  used  to  generate  the  Ni-rich 
part  of  the  Al-Ni  phase  diagram  based  as  shown  in 
Fig.  10.  Here  the  two-phase  region  is  narrowed  signifi¬ 
cantly,  with  the  y'/ (y'  +  y)  boundary  moving  to  higher  Ni 
concentrations  and  the  high-temperature  portion  of  the 
(y'  -h  y)/y  boundary  moving  to  slightly  lower  Ni  concentra¬ 
tions  relative  to  the  first  cluster  expansion.  Evaluating  the 
low-temperature  sharp  interfacial  energy  using  this  cluster 
expansion  produces  energies  of  42,  46  and  39  mJ  for 
the  (100),  (110),  and  (111)  boundaries,  respectively. 
This  is  in  much  better  agreement  with  energies 
calculated  directly  from  VASP  (Section  3.1),  41,  55  and 
48  mJ  m“^  for  the  (100),  (110)  and  (111)  boundaries, 
respectively. 

Fig.  11  shows  the  results  of  integrating  interfacial  excess 
energy  (Eq.  (1))  produced  with  CE2  for  the  three  low-index 
boundaries  in  order  to  produce  the  lEFE.  As  shown  in  the 
figure  these  energies  drop  rapidly  with  increasing  tempera¬ 
ture,  becoming  constant  and  fall  in  a  narrow  energy  range 
(9-12  mJ  m“^)  for  temperatures  above  600  K.  The  second 
cluster  expansion  produces  an  average  lEFE  of  approxi¬ 
mately  1 1  mJ  m~^  for  the  three  low-index  boundaries, 
which  is  very  good  agreement  with  LSW  analysis  published 
in  the  mid-1990s  [8]. 

The  two  cluster  expansions  exhibit  a  rapid  decrease  in 
the  IEEE  with  temperature,  with  CE2  producing  a  narrow 
window  of  the  three  free  energies  over  most  of  the  temper¬ 
ature  distribution.  At  high  temperatures,  above  1000  K,  the 
two  approximations  give  similar  results  but  with  different 
slopes  for  the  (1 1 0)  and  (110)  lEFEs.  Because  of  the  slow 


Atomic  Fraction  Ni 

Fig.  10.  The  Ni-rich  phase  diagram  produced  by  the  second  cluster 
expansion  (CE2).  Solid  lines  are  from  CALPHAD  calculations  using  the 
PANDAT  method  [30,40]  and  the  dashed  lines  are  cluster  expansion 
predictions  using  Monte  Carlo  methods.  By  initializing  the  ground  state 
search  within  AT  AT  using  fee  Ni,  LI2  and  the  relaxed  (1 00),  (1 1 0) 

and  (111)  IFBs  the  y'/{y-\-y')  phase  boundary  is  shifted  to  higher  Ni 
compositions  relative  to  CEl.  The  {y  +  y')/y  phase  boundary  remains 
relatively  unchanged  compared  to  CEl.  The  vertical  dotted  line  in  the 
figure  marks  the  Ni-15A1  (at.%)  composition. 


Temperature  (K) 

Fig.  11.  The  lEFE  as  a  function  of  temperature  for  the  (1 00),  (1 1 0)  and 
(111)  y-y'  boundaries  calculated  using  CE2.  The  lEFE  is  derived  by 
numerically  integrating  Eq.  (1)  using  the  piecewise  continuous  polynomial 
fit  to  the  lEFE.  The  solid  and  dashed  lines  are  the  predicted  lEFEs,  and 
are  in  good  agreement  with  lEFEs  derived  from  earlier  coarsening 
experiments  of  Ardell.  Open  and  filled  circles  are  from  coarsening  data  fit 
to  LSW  theory  [8]  and  the  filled  square  is  from  fits  using  trans-interface- 
diffusion-controlled  theory  [10]. 


kinetics  of  these  alloy  systems  micro  structural  evolution 
simulations  would  typically  be  in  the  range  of  800-1450  K. 

Given  the  key  role  the  y'  particles  play  in  the  high-tem- 
perature  strength  of  Ni-based  superalloys  it  is  worth  spec¬ 
ulating  on  mechanisms  that  could  stabilize  the  precipitates 
to  higher  temperatures.  As  the  temperature  rises  the  con¬ 
figurational  entropy  in  Ni  increases  which  acts  to  reduce 
the  lEFE.  Given  the  strong  evidence  of  the  role  of  the 
local  magnetic  moment  of  Ni  in  increasing  the  IFB  energy 
it  is  reasonable  to  assume  that  introducing  atomic  species 
with  strong  local  magnetic  moments,  that  preferentially 
favor  the  y  phase,  could  stabilize  the  IFB  to  higher 
temperatures. 

5.  Summary 

First-principles  methods  were  used  to  assess  the  struc¬ 
ture,  energies  and  free  energies  of  low-index  y-y'  bound¬ 
aries  in  the  Ni-Ni3Al  system.  We  verify  the  previous 
work  of  Price  and  Cooper  that  within  the  GGA  approxi¬ 
mation  these  IFB  energies  are  positive  only  when  one 
includes  the  spin  degrees  of  freedom.  The  relative  energies 
of  the  low-temperature  chemically  sharp  IFBs  reflect  the 
number  of  Ni  spin  moments  that  are  suppressed  by  the 
introduction  of  the  interface.  This  effect  appears  to  be  rel¬ 
atively  short  range  for  the  (100)  and  (110)  IFBs  and  longer 
range  for  the  (111)  boundary.  However,  to  a  first-order 
approximation  the  largest  contribution  to  the  energy  is 
due  to  the  spin  suppression  in  the  first  all-Ni  layer  adjacent 
to  the  boundary. 

Spin-polarized  Ni-Al  reference  configurations  are  also 
needed  when  developing  effective  cluster  interaction  coeffi¬ 
cients  for  use  in  a  cluster  expansion  description  of  the  Ni- 
rich  portion  of  the  binary  phase  diagram.  Cluster  expan¬ 
sions  based  on  spin-averaged  calculations  produce  at  least 
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one  unexpected  phase  (NijAl)  in  the  Ni-rich  portion  of  the 
binary  phase  diagram  that  is  stable  at  high  temperatures. 
This  result  is  consistent  with  the  previous  density  func¬ 
tional  calculations  by  Wolverton  and  Zunger  [12],  who 
found  that  magnetization  destabilizes  the  NijAl  phase  rela¬ 
tive  to  Ni  and  Ni^Al. 

Using  a  simple  cluster  expansion  in  conjunction  with 
Monte  Carlo  methods  we  find  the  widths  of  the  low-index 
IFBs  in  the  range  of  1. 0-3.0  nm  depending  on  the  bound¬ 
ary  and  the  temperature.  The  width  of  the  (111)  IFB  is 
generally  ^  l.l  times  larger  than  the  (100)  and  (110) 
boundaries  at  any  given  temperature.  The  rate  of  increase 
in  width  as  a  function  of  temperature  ranges  from  0.2  to 
0.3  nm  per  100  K.  Given  the  relatively  narrow  IFB  widths 
we  conclude  that  phase  evolution  methods  that  explicitly 
model  the  boundaries  will  be  limited  to  relatively  small 
simulation  volumes.  This  approach  also  predicts  a  small, 
temperature-dependent,  shift  in  the  centers  of  the  bound¬ 
aries  as  defined  through  the  concentration  and  LRO 
parameter.  The  shift  is  smaller  than  that  observed  in 
Rene  88DT  using  high-resolution  TEM  imaging  methods 
[32]. 

The  lEFEs  are  also  strongly  temperature  dependent, 
falling  to  10  mJ  at  1000  K.  The  largest  entropic  con¬ 
tributions  to  the  free  energy  are  from  the  configurational 
degrees  of  freedom  while  the  vibrational  and  electronic 
terms  are  smaller  for  the  disordered  IFBs.  The  lEFE  has 
been  studied  using  two  different  cluster  expansions  and 
the  general  trends  are  insensitive  to  the  specifics  of  ECIs. 
The  results  for  the  binary  Ni-Al  system  are  encouraging. 
Given  the  chemical  complexity  of  commercial  Ni-based 
superalloys  studies  of  the  effects  of  high  concentrations  of 
ternary  elements  (e.g.  Cr  or  Co)  on  the  IFBs  would  be 
highly  desirable. 
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